Consistent LDA'+DMFT approach to electronic structure of transition metal oxides: 
charge transfer insulators and correlated metals. 
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We discuss the recently proposed LDA'+DMFT approach providing consistent parameter free 
treatment of the so called double counting problem arising within the LDA+DMFT hybrid com- 
putational method for realistic strongly correlated materials. In this approach the local exchange- 
correlation portion of electron-electron interaction is excluded from self consistent LDA calculations 
for strongly correlated electronic shells, e.g. d-states of transition metal compounds. Then the 
corresponding double counting term in LDA'+DMFT Hamiltonian is consistently set in the local 
Hartree (fully localized limit - FLL) form of the Hubbard model interaction term. We present the 
results of extensive LDA'+DMFT calculations of densities of states, spectral densities and opti- 
cal conductivity for most typical representatives of two wide classes of strongly correlated systems 
in paramagnetic phase: charge transfer insulators (MnO, CoO and NiO) and strongly correlated 
metals (SrV0 3 and Sr 2 Ru0 4 ). It is shown that for NiO and CoO systems LDA'+DMFT quali- 
tatively improves the conventional LDA+DMFT results with FLL type of double counting, where 
CoO and NiO were obtained to be metals. We also include in our calculations transition metal 
4s-states located near the Fermi level missed in previous LDA+DMFT studies of these monooxides. 
General agreement with optical and X-ray experiments is obtained. For strongly correlated metals 
LDA'+DMFT results agree well with earlier LDA+DMFT calculations and existing experiments. 
However, in general LDA'+DMFT results give better quantitative agreement with experimental 
data for band gap sizes and oxygen states positions, as compared to the conventional LDA+DMFT. 

PACS numbers: 71.20.-b, 71.27+a, 71.28+d, 74.25.Jb 



I. INTRODUCTION 

During last decade the LDA+DMFT method (local 
density approximation + dynamical mean-field theory) 
became probably the most powerful tool to calculate elec- 
tronic structure of real strongly correlated materials [l|- 
[||. Typically this approach consists of two computation 
steps. First, LDA calculations are exploited to obtain 
the non-interacting Hamiltonian H LBA which describes, 
rather accurately, the kinetic energy (and to some extent 
takes into account electronic interactions). At the second 
step the local Coulomb (Hubbard) interaction iJ Hub is 
introduced into the lattice problem defined by H LBA for 
those electronic shells which are supposed to be strongly 
correlated. Thus obtained generalized Hubbard model is 
solved numerically using DMFT. Some attempts to or- 
ganize a feedback from DMFT step back to LDA calcu- 
lations to achieve fully self-consistent LDA+DMFT are 
also known and may be important for some physical prob- 
lems 0- 

The double counting problem arises in the standard 
LDA+DMFT, because some portion of local electron- 
electron interaction for correlated shells is actually ac- 
counted for within i? LDA . To avoid this double counting 
it is necessary to subtract a certain correction term H DC 
from H LDA . Then the formal LDA+DMFT Hamiltonian 
is written as: 
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In orbital space H DC is the diagonal matrix with non 
zero and equal matrix elements Ed c for these atomic 
shells assumed to be strongly (e.g. d or / shells or their 
subshclls). This becomes more transparent if we con- 
sider the corresponding Green's function for the Hubbard 
model: 

GijQtE) = [(B-/i)7-^. DA (k)-(S(kS)-S dc )^ d ]- 1 , 

(2) 

where I is the unity matrix in the orbital space, is the 
chemical potential and S(k£J) is the self-energy corre- 
sponding to local Coulomb (Hubbard) interaction, 
denotes matrix inversion, while index d denotes corre- 
lated states for which Coulomb (Hubbard) interaction is 
taken into account. 

From Eq. ([2]) one can see that in case of iJ LDA con- 
taining only the contribution of interacting d - orbitals, 
Edc reduces to trivial renormalization of the chemical 
potential fj.. Then, strictly speaking, there is no double 
counting problem at all. Because of this many of early 
works (listed e.g. in reviews [2 H-Q , except probably 
the first paper on LDA+DMFT 1J and few others) just 
dropped the double counting correction term. Only after 
the LDA+DMFT community started the active studies of 
multiband if LDA Hamiltonians with both correlated and 
non correlated states included, the problem of correct 
implementation of H DC became important. Now there 
are dozens of works devoted to multiband LDA+DMFT 
studies. Important classes of materials investigated can 
be listed as: 
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1. Transition metal oxides (LaTi03, (Sr,Ca)V03, 
V 2 3 , V0 2 , Cr0 2 , LaMn0 3 , NiO, MnO, CoO, 
FeO, LaCo0 3 , TiOCl, Tl 2 Mn 2 7 , LaNi0 3 , 
(Ca,Sr) 2 Ru0 4 , Na . 3 CoO 2 ); 

2. Elemental transition metals and non-oxide transi- 
tion metal compounds (Cr, Mn, Fe, Ni, Co, mul- 
tilayers (CrAs)/(GaAs), NiMnSb, Co 2 MnSi, CrAs, 
VAs, ErAs, Ni(S,Se) 2 , KCuF 3 ); 

3. Elemental /-electron materials and their com- 
pounds (Ce, Pu, Am, Ce 2 3 , Pu 2 3 , USe, 
UTc, PuSe, PuTe, PuCoGa 5 , URu 2 Si 2 , Celrln 5 , 
CeCoIn 5 , CeRhIn 5 ); 

4. Nano materials (Ni-Cu nano contacts and nano 
electrodes); 



JT]) is given by: 



5. High temperature copper 
((Sr,La) 2 Cu0 4 , (Pr,Ce) 2 Cu0 4 , 
etc.); 



superconductors 
Bi 2 Ca 2 SrCu0 8 



6. Superconducting iron pnictides (LaFeAsO, Ce- 
FeAsP, LiFeAs, BaFe 2 As 2 , etc.). 

These systems show a large variety of physical effects. 
Among them there are strongly correlated metals, Mott 
and charge transfer insulators, ferromagnets and antifcr- 
romagnets, superconductors, etc. However, up to now 
there is no universal and unambiguous expression for 
H DC , and different formulations are used for different 
classes of materials. 

In this paper we present the results of extensive ap- 
plication of our recently proposed LDA'+DMFT [ljj ap- 
proach to charge transfer insulators MnO, CoO and NiO 
and strongly correlated metals SrV03 and Sr 2 Ru0 4 , con- 
fronted to conventional LDA+DMFT results and some 
experiments. The manuscript has following structure. 
In Sec. II we present an overview of different defini- 
tions of the H DC . The novel consistent LDA'+DMFT 
method is described in Sec. III. LDA and LDA' band 
structures, total and partial densities of states, spectral 
density maps and optical conductivity LDA'+DMFT re- 
sults for prototype charge transfer insulators MnO, NiO 
and CoO are presented in Sec. IV and compared with the 
results of conventional LDA+DMFT. These results are 
further compared with experimental data on X-ray spec- 
troscopy and optical conductivity. In Sec. V we discuss 
LDA and LDA' band structures for correlated metallic 
systems prototypes SrV03 and Sr 2 Ru0 4 are presented. 
Then LDA+DMFT and LDA'+DMFT results arc com- 
pared with each other and with experimental photoemis- 
sion and absorption spectra. Finally we end up with the 
Conclusion (Sec. VI). 



II. REVIEW OF DIFFERENT FORMULATIONS 



FOR H 



DC 



To derive an expression for H DC let us examine iJ LDA 
and H Hub terms Eq. (TTJ . LDA part of the Hamiltonian 
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(3) 



where A is the Laplace operator, m e the electron mass, 
e the electron charge, and 



(4) 
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denote the one-particle potential due to all ions i with 
charges eZi at given positions Ri, and the electron- 
electron interaction, respectively. 

The £"x C DA (p(r)) in © is some function of local charge 
density, which approximates the true exchange correla- 
tion functional E^ c [p] of density functional theory within 
local density approximation Q. The explicit expression 
for El? A (p(r)) is usually derived from perturbation the- 
ory Q or numerical simulations fiol ] of the "jellium" 
model with Vi on (r) = const. To obtain the value of lo- 
cal charge density one should choose some basis set of 
one-particle wave functions <fi (e.g. to do practical cal- 
culations and explicitly express matrix elements of the 
Hamiltonian ((SJ), so that p(r) is written as: 



N 



p(r) = 5> 4 (r)| 2 



(5) 



Hubbard-like (local) interaction term including direct 
Coulomb interaction and exchange Coulomb interaction 
contributions in the density-density form is written as: 



m i 

+ E E E^'-^' J ) 



Hma ' L im' a' 



(6) 



Here, the index i enumerates lattice sites, m denotes or- 
bitals, and a the spin. The U represents local intra- 
orbital Coulomb repulsion and J - z-component of 
Hund's rule coupling between the strongly correlated 
electrons (e.g. d-states, enumerated by i = and I = Id)- 
Rotational invariance then fixes the local inter-orbital 
Coulomb repulsion U' — U — 2J 11 . The values of U 
and J are obtained usually from constrained LDA pro- 
cedure [lil ]. One can get numerically exact solution of 
the Hubbard Hamiltonian (simplified kinetic term plus 
jjHub t erm ) within DMFT approximation. 

Hamiltonian -£?lda contains local electron-electron cor- 
relations through the exchange correlation energy (taken 
in the form valid for uniform electronic gas) and density- 
density contribution of the Hartree term. In its turn, 
DMFT provides the numerical solution of the Hubbard 
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model (exact in infinite dimensions). Thus it is clear 
that before plugging -Hlda into DMFT lattice problem 
([2]), one must subtract certain double counting correc- 
tion term H DC from Hlda- The double counting prob- 
lem arises because there is no explicit microscopic or di- 
agrammatic relation between the model (Hubbard like) 
Hamiltonian approach and LDA. There is apparently no 
possibility to give a rigorous expression for H in terms 
of U, J and p. Thus, several ad hoc expressions for H DC 
and approaches to treat the double counting problem ex- 
ist in the current literature. Below we briefly discuss 
some of these derivations. 

Perhaps for the first time problem of double count- 
ing appeared within an attempt to merge LDA and the 
Hubbard model within the LDA+U method [HJ], where 
was initially postulated the so called "around mean-field" 
(AMF) definition of H DC . This definition comes from 
an assumption that LDA is a kind of "mean-field" solu- 
tion of the Hubbard-like problem Eq. Later on the 
definition of Ref. was generalized for spin dependent 
(LSDA) case (and even more general - with matrix form 
of Coulomb interaction). After this spin dependent gen- 
eralization corresponding AMF expression can be given 
as: 

h amf = ^U^2n da (n d - -J^2n da (n dc - n°) 

(7) 

with the average occupancies n° = 2 (2i+i) Em a n ma: 
n a = (2i+i) Em n ™a and total number of electrons 
on interacting orbitals (per spin projection) n da — 

J2m n Udm<r = Em(™''d™) and n d = E CT n <fo- originally 
supposed to be found from LDA calculations. The draw- 
back of AMF is the equal occupancy of all orbitals which 
is not correct even for weakly correlated systems because 
e.g. of crystal field splitting. However, a couple of the 
modern LDA+DMFT works reported the reasonable re- 
sults with AMF-like double counting correction term. 
Apparently, the AMF double counting correction works 
rather well for moderately correlated metallic systems. 
Some modifications of ([7J were given in Refs. [16| and 
applied to LDA+DMFT calculations for charge transfer 
insulators. 

Later on the fully localized (or atomic) limit (FLL) 
expression for H 6 was introduced in Refs. 0171 (with 
first application to LDA+DMFT calculations in Ref. pj): 



H^L = fJn d {n d 



n da (n da - 1). (8) 



The Eq. © actually represents the Hartree decoupling of 
the Hubbard model interaction term (|6]) — decoupling of 
the density-density term hihj and not full four operator 
term ctcjc c;. Thus strictly speaking in Eq. ([SJ there 
is no Fock type of contribution since Hund exchange is 
presented in Eq. ([6]) in the density-density form, although 
Hund coupling value J has "exchange nature" . Quite 



often it is misinterpreted as due to the "true" Hartree- 
Fock decoupling of clcjC Ci term. 

The FLL expression in the context of LDA+DMFT 
calculations was used in the majority of modern works. 
It works reasonably good for both metallic and insulat- 
ing strongly correlated materials. Recent ly som e mod- 
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ifications of FLL were proposed in Refs. 
ically these modifications are used to quantitative im- 
provements of LDA+DMFT results for particular com- 
pounds. Some kind of AMF and FLL "hybrid scheme" 
was used in Ref. [2(J for a — Fe. 

Alternative way to derive or guess the H DC term is to 
express it through the characteristics of intrinsic single 
DMFT impurity problem, such as impurity self-energy 
^mm' or impurity Green's function G 1 ™^, . A popular 
way is to define double counting energy as a static part 
of the impurity self-energy [2l| : 



E dc = iT 1V (Sr P (0)). 



(9) 



Some of LDA+DMFT papers used this definition in cal- 
culations of metallic magnetic and non-magnetic systems. 
From the very beginning this type of double counting 
correction was also exploited within the GW+DMFT ap- 
proach [22j . 

Hartree energy can be determined from LDA+DMFT 
self-energy as its real part in the high frequency limit 
value. In Ref. [25| it was proposed to use thus defined 
Hartree energy as a double counting correction, using the 
constraint 



ReTr(E^,(^))=0, 



(10) 



where lo^ is the highest Matsubara frequency (used in 
calculations). Physically similar definition of double 
counting term E dc = E(w — > oo) was successfully applied 
to metallic ferromagnetic SrCoQ^p3l|. 

For metallic systems it was suggested to fix the double 
counting correction by equating the number of particles 
of non-interacting problem and impurity problem as ex- 
pressed via corresponding Green's function [24| : 



Tr G l 



TrG^(/3), 



(11) 



where G^ m i is local non interacting Green function. 
Some of LDA+DMFT works treated double counting en- 
ergy E dc as a free parameter. The authors of Ref. [25[ 
found that most of described H DC terms proposed in the 
literature are not completely satisfactory in the case of 
charge transfer insulator NiO and proposed a numerical 
way to define the necessary double counting correction. 

Another possible solution of the double counting 
problem is to perform Hartree+DMFT or Hartree- 
Fock+DMFT calculations |26(. While performing 
Hartree-Fock band structure calculations for real ma- 
terials we do know exactly what portion of interaction 
is included. Since diagrammatic expression for Hartree 
or Hartree-Fock terms are well known, one can calcu- 
late them directly and get double counting correction 
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energy explicitly. However, up to now we are unaware 
of any Hartree+DMFT or Hartree-Fock+DMFT calcu- 
lations for real materials. 

Completely independent branch of ab initio DMFT cal- 
culations is GW+DMFT method, which uses instead of 
density functional theory the so called chain of Hedin 
equations truncated in a simplest manner by the neglect 
of vertex corrections (for review see Ref. [22.27)]. Be- 
cause of purely diagrammatic nature of GW there is a 
natural way to calculate the local part of corresponding 
Hartree contribution, which can be used as the double 
counting correction term for GW+DMFT. [27|. 



III. CONSISTENT LDA'+DMFT APPROACH 

Recently we proposed the LDA'+DMFT approach, 
which defines consistent parameter free way to avoid the 
double counting problem [l2j ]. The main idea is to ex- 
clude explicitly exchange-correlation energy from self- 
consistent LDA calculations only for correlated bands. 
As described above main obstacle to express double 
counting term exactly is exchange correlation E^ A (p(r)) 
portion of interaction within LDA. So it seems somehow 
inconsistent to use it to describe correlation effects in nar- 
row (strongly correlated) bands from the very beginning, 
as these should be treated via more elaborate schemes 
like DMFT. To overcome this difficulty for these states, 
we propose to redefine charge density (J5)) in E^ A as 
follows: 



P'(r) 



(12) 



excluding the contribution of the density of strongly cor- 
related electrons. 

In principle E^ A is not an additive function of charge 
density. Thus splitting of charge density into two parts 
may lead to some loss of hybridization between correlated 
and uncorrelated states. However as we show below this 
approximation is rather good. Later on we will see that 
LDA' bands practically do not change their shape with 
respect to LDA ones for all considered systems. That 
tells us that "hybridization" is almost not affected by 
LDA'. The main effect is increase of splitting between 
oxygen 2p and metal 3d states. It comes from more re- 
pulsive potential appearing in the LDA' case since part 
of exchange correlation energy is excluded there. 

Then this redefined p'(r) (JT^J) is used to obtain ££ C DA 
and perform the self-consistent LDA' band structure cal- 
culations for correlated bands. This procedure leaves out 
of interaction for correlated states on the LDA' stage 
just the Hartree contribution ([3]). Thus, double counting 
correction term should be consistently taken in the form 
of the Hartree like term, given by Eq. ©. This Hp£ L 
definition also does not have any free parameters. Actu- 
ally, our approach is in precise correspondence with the 
standard definition of correlations, as interaction correc- 
tions "above" the Hartrcc-Fock. At the same time all 



other states (not counted as strongly correlated) are to 
be treated with the full power of DFT/LDA and full p in 

rnLDA 
fc xc ■ 

Despite the fact that the LDA'+DMFT method is ap- 
parently most consistent with the use of FLL type of dou- 
ble counting, in principle all mentioned above definitions 
oiH DC can also be exploited within LDA'+DMFT. Also 
there is another "degree of freedom left" - the occupancy 
ndi used in FLL equation, can be obtained either from 
LDA or LDA' results, or it can be calculated self consis- 
tently during the DMFT loop. We used all these variants 
in our calculations for different compounds presented be- 
low. Corresponding values of Ed c listed in Table HI No- 
tations arc: FLL(SC) for self consistently calculated rid 
and for n d from LDA or LDA' - FLL (LDA). In general 
FLL(SC) and FLL(LDA) do not differ very much from 
each other, except for the case of CoO (see below). How- 
ever, FLL(SC) gives slightly better agreement with ex- 
periments. Most Figures presented below arc plotted for 
the FLL(SC) case. We observed that FLL(SC) calcula- 
tions require more computational time than FLL(LDA). 

Thereby our consistent LDA'+DMFT approach is 
a kind of compromise between Hartree-Fock and 
DFT/LDA starting points to be followed byDMFT cal- 
culations. It was demonstrated in Ref. [l2| that this 
LDA'+DMFT method works perfectly for insulating NiO 
system, directly producing charge transfer insulator solu- 
tion, while conventional LDA+DMFT (with FLL) gives 
metallic solution (cf. Ref. plj). 



IV. CHARGE TRANSFER INSULATORS 

A. LDA and LDA' band structures 

Typical examples of charge transfer insulator (CTI) 
materials are transition metal monoxides MnO, CoO and 
NiO. These oxides have rock salt crystal structure with 
lattice parameters a=4.426A, 4.2615A and 4.1768A cor- 
respondingly. To obtain LDA and LDA' band structures 
for MnO, CoO and NiO the linearized muffin-tin orbitals 
(LMTO) basis set [HI was used. In the corresponding 
program package TB-LMTO v.47 the E^ A was taken 
in von Barth-Hedin form [§]. Total and partial densities 
of states (DOS) together with band dispersions can be 
seen in Fig. [1] for LDA (dashed lines) and LDA' (solid 
lines). From top to bottom on Fig. [T] there are MnO, 
CoO and NiO systems. As reported earlier for NiO [l2[ 
LDA' approach changes charge transfer energy \Ed — E p \, 
where E d and E p are, roughly speaking, one electron en- 
ergy positions of transition metal 3d and 0-2p bands. In 
Fig. [1] the same tendency for MnO and CoO oxides is 
seen. For MnO it increases about 0.5 eV and for CoO 
about 1 eV, similar to NiO. Almost rigid shift of 0-2p 
bands down in energy is observed here, while transition 
metal 3d states remain almost the same near the Fermi 
level. 

One should mention that (to our knowledge) transi- 
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FIG. 1: LDA (dashed lines) and LDA' (solid lines) densities 
of states (DOS) and band dispersions for MnO (upper row), 
CoO (middle row) and NiO (lower row). Fermi level is zero. 



tion metal 4s states were never included previously into 
LDA+DMFT calculations for these transition metal ox- 
ides. Apparently, this happened because they were rea- 
sonably assumed to be weakly correlated and thus pro- 
jected out from corresponding LDA Hamiltonian. How- 
ever, transition metal 4s states are rather close to the 
Fermi level for LDA bands and getting even closer for 
LDA' ones. They can be seen on Fig. [T] as lowest unoc- 



cupied states which are touching the Fermi level for MnO 
near T-point and less than leV above the Fermi level for 
CoO and NiO. 



B. LDA+DMFT and LDA +DMFT spectral 
functions 



Everywhere in this paper we employ Hirsh-Fye quan- 
tum Monte- Carlo algorithm [29j as impurity solver for 
DMFT equations. To set up DMFT lattice problem we 
use corresponding LDA and LDA' Hamiltonians, which 
include all states (without any projecting, as was done 
e.g. Ref. [13 )■ Inverse temperature was taken fj = 
5cV _1 , with 80 time slices for NiO, while for MnO and 
CoO we used (3 = 10cV _1 with 120 and 160 time slices 
respectively. Monte Carlo sampling was done with 10 6 
sweeps. The use of rather high temperatures does not 
lead to any qualitative effects in the results, allowing 
us to avoid unnecessary computational efforts. Parame- 
ters of Coulomb interaction were chosen as typical for 
MnO, CoO and NiO [THi]: ^=8 cV and J=l cV. 
Both FLL(SC) and FLL(LDA) double counting defini- 
tions were applied for all materials. Respective E ( i c val- 
ues are given in Tabic HI 

To obtain DMFT(QMC) densities of states (DOS) at 
real energies, we employed the maximum entropy method 
(MEM) [3(|. Then one can get DMFT self-energy on the 
real frequency axis by using Pade approximants for an- 
alytical continuation. Further on it was checked that 
"Pade" DOS'es are identical to "MEM" DOS'es. Once 
£(w) is obtained, one can input it into Eq. ([2]) and obtain 
the spectral density function A(k,ui) = — -ImG(k, ui). 
Corresponding maps of spectral density functions, rep- 
resenting effective band structure of these compounds, 
are given in Fig. [2] Left column of Fig. [5] presents 
LDA+DMFT results and the right one - LDA'+DMFT 
for MnO (upper panels), CoO (middle panels) and NiO 
(lower panels). 



C. LDA+DMFT and LDA'+DMFT DOS 

In Fig. |3] we present densities of states obtained by 
LDA+DMFT (dashed lines) and LDA'+DMFT (solid 
lines). The left panel corresponds to MnO, middle one 
to CoO and the left one to NiO. Upper row shows to- 
tal densities of states, while in other rows we show the 
contributions of the most important electron states - t2 g 
and e g subshells for 3d transition metal, oxygen 2p states 
and also transition metal 4s states. 

First we focus on MnO case which is perhaps the sim- 
plest one among these three. The 0-2p states are located 
between -9 and -4 eV (see Figs. and [3]). Then comes 
lower Hubbard band, which consists of Mn-3d ti g and e g 
contributions at -4eV and -2.3 eV correspondingly. On 
spectral function maps LHB is rather wide non disper- 
sive band at these energies. Then we see the so called 
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FIG. 2: (Colour online) Comparison of LDA+DMFT (left column) and LDA'+DMFT (right column) calculated spectral density 
functions for MnO (upper row), CoO (middle row) and NiO (lower row), with FLL(SC) double counting correction. Fermi level 
is zero. 



Zhang-Race band - the bound state which appears when 
strongly interacting band is hybridized with charge reser- 
voir. This band can be seen as a peak at -1.5 eV in 0-2p 
states together with Mn-3d e g states. Then, between the 
Zhang- Race band and the upper Hubbard band there is a 
gap for Mn-3d states about 3.5 cV for both LDA+DMFT 
and LDA'+DMFT, which agrees pretty well with exper- 
imental spectra (see below). UHB is located above 4 eV 
and, where ti g and e g contributions can not be separated 
in energy. 

Spectral density map of Fig. [2] (upper row) show some 
rather well defined band of MnO, which touches the 
Fermi level in the T-point. This band is nothing else 
but Mn-4s. It is seen from Fig.[3j that most of the Mn-4s 
spectral weight is actually well above 5 eV. Below there 
is some rather low intensity tail, which goes through the 



gap between the upper Hubbard band and the Zhang- 
Rice band. Its intensity is at least one order of magnitude 
lower, than intensities of other contributions to DOS. 

Consider next CoO (middle row of Fig. [2] and mid- 
dle panel of Fig. [3]). We see that both LDA+DMFT 
and LDA'+DMFT results are quite similar. There is 
some difference in the UHB, where Co-3d ti g and e g con- 
tributions can now be separated and in Fig. [3] two al- 
most nondispersive bands around 2 and 3 eV above the 
Fermi level are clearly seen. The gap between Zhang- 
Rice band and UHB is about 0.5 eV larger (about 4 eV) 
for LDA'+DMFT results. 

One should note, that LDA+DMFT calculation with 
FLL(LDA) double counting produces the metallic so- 
lution for CoO, as seen from Fig. 21 which qualita- 
tively contradicts the experiments. On the contrary, 
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FIG. 3: Comparison of LDA+DMFT (dashed lines) and LDA'+DMFT (solid lines) densities of states for MnO (left panel), 
CoO (middle panel) and NiO (right panel), with FLL(SC) double counting correction. Fermi level is zero. 
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FIG. 4: Comparison of LDA+DMFT (dashed lines) and 
LDA'+DMFT (solid lines) calculated densities of states for 
CoO with FLL(LDA) double counting correction. Fermi level 
is zero. 



LDA'+DMFT gives the correct insulating state. 

Note, that both in CoO and NiO the behavior of 4s 
bands is similar to that discussed above for the case of 
MnO. Spectral density maps of Fig. [2] show the presence 
of these bands within the charge transfer gap, though 
the partial density of states due to these bands within 
the gap is almost negligible (cf. Fig. [3]) . 

To sum up, we stress that both MnO and CoO 
within LDA'+DMFT arc consistently demonstrated to 
be charge transfer insulators (in contrast to the conven- 
tional LDA+DMFT in the case of CoO). The similar be- 
havior was obtained earlier for NiO in Ref. [l2[ . Here we 
presented more complete LDA'+DMFT results for NiO, 
with both FLL(LDA) and FLL(SC) double counting cor- 
rections. Conventional LDA+DMFT calculations pre- 
dict NiO to be metallic in contrast to experiment, while 
LDA'+DMFT gives charge transfer insulating solution 
for NiO for both FLL(LDA) and FLL(SC) double count- 
ing correction. All other features of NiO LDA'+DMFT 
band structure are quite similar to MnO and CoO com- 
pounds described above. 



D. LDA+DMFT and LDA'+DMFT optical 
conductivity 

Metallic or insulating behavior can be explicitly 
demonstrated by calculations of optical conductivity. Be- 
low we present our results for optical conductivity be- 
havior of MnO, CoO and NiO within LDA+DMFT and 




FIG. 5: Comparison of experimental (circles, stars) and calculated LDA+DMFT (dashed lines) and LDA'+DMFT (solid lines) 
optical conductivities for MnO (left panel), CoO (middle panel) and NiO (right panel). 



LDA'+DMFT, allowing us also to analyze the influence 
of transition metal 4s states on dielectric properties of 
these oxides. In our calculations we used the following 
expression for optical conductivity, valid in DMFT [3~j| : 

oo 

^e 2 f , /(e) - f(e - u) 

cr xx {oj) = — / de^ — ^ '- x 

Zna J oj 

oo 

* gW<.Mf<«-«>(u) 

ijka \ / \ / 

Here e is electron charge, a is the lattice constant of 
corresponding compound, /(e) - Fermi function, e^ - 
band dispersion, A l 2 (e) corresponding (LDA+DMFT or 
LDA'+DMFT) spectral density function matrix are 
the band indices). During our calculations we found that 
main contribution to optical conductivity is due to intra- 
orbital optical transitions. Inter-orbital optical transi- 
tions give less than 5% of optical conductivity intensity 
in frequency range used in our calculations. Also in the 
present work we neglect possible effects due to optical 
matrix elements. Calculated theoretical curves ob- 
tained in conventional LDA+DMFT (dashed line) and 
within LDA'+DMFT (solid line) are presented in Fig. [5] 
for MnO (left panel), CoO (middle panel) and NiO (right 
panel) . 

From Fig. E] we see, that within LDA'+DMFT (solid 
line) all materials are insulators. Despite the presence 
of transition metal 4s states close to the Fermi level, 
possible Drude peak due to these states is not observed. 
Conventional LDA+DMFT optical conductivity for NiO 
shows typical metallic behavior, as discussed earlier in 
the context of DOS behavior. 

Now we compare our theoretical results with available 
experimental data (with an exception of MnO, where we 
are not aware of any experimental results) [32||33| . In 
Ref. [lij only experimental data for optical constants 
n{uS) and k(cj) were presented. The optical conductiv- 
ity in units of j- (which is about 5.8x 10 3 il _1 cm _1 for 
given monooxides) can be recalculated from these data 

using as cr(ui) = %!l woT 1 ^, where a is fine struc- 
ture constant, a - lattice constant and c - speed of light. 



Corresponding curves are shown in Fig. [5] by stars. For 
NiO there are more recent experimental data of Ref. [33| , 
shown with circles. One observes that below the leading 
absorption edge for CoO and NiO there exist rather long 
absorption tails with low intensity. We associate these 
tails with contribution of Co and Ni 4s states. For NiO 
the overall agreement of LDA'+DMFT results with ex- 
perimental data is quite satisfactory. For CoO theoretical 
absorption edge is about 1 eV lower than experimental 
one. However, this can probably be corrected introduc- 
ing the larger value of Coulomb interaction U. Recent 
constrained RPA study produced it to be 10.8 eV [23| . 
in contrast to 8 eV used in our calculations. 



E. Comparison of LDA+DMFT and LDA'+DMFT 
results with X-ray experiments 

Now we compare our results for DOS with XPS and 
BIS experiments of Refs. [Mil. In Fig. © LDA+DMFT 
(dashed lines) and LDA'+DMFT (solid lines) valence 
and conduction bands spectra are directly compared with 
spectra for MnO (upper panel), CoO (middle panel) and 
NiO (lower panel) . Theoretical spectra were obtained by 
multiplication of DOS by Fermi distribution and Gaus- 
sian broadening with experimental temperature and res- 
olution. 

General structure of spectra is similar for all three com- 
pounds. From -14 to -4 eV there are 0-2p states, then 
comes lower Hubbard band at about -3 eV. On the high 
energy slope of the LHB there we can see a shoulder-like 
structure, which is nothing else but Zhang-Rice band. 
Around the Fermi level there is insulating gap. The 
size of the gap is very well reproduced for MnO by both 
LDA+DMFT and LDA'+DMFT. For CoO it looks like U 
value chosen is a bit too small (as discussed earlier) , how- 
ever LDA'+DMFT spectra gives gap size closer to the 
experiment. For NiO conventional LDA+DMFT gives 
metallic solution, while LDA'+DMFT produces CTI so- 
lution with correct energy gap size. Experimental po- 
sitions of the upper Hubbard bands arc rather well de- 
scribed by LDA'+DMFT. Since experimental data for 
NiO goes far above the Fermi level one can identify these 
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On Fig. [6] one can see that experimental conduction 
band low energy threshold has a rather long low intensity 
tail which goes down to the Fermi level. Because of that 
there is some asymmetry of the gap. We suggest that this 
asymmetry of the gap originates from transition metal 
4s states, which touch the Fermi level from above, as 
described earlier. 



V. STRONGLY CORRELATED METALS 



LDA and LDA' band structure 
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FIG. 6: Comparison of LDA+DMFT (dashed lines) and 
LDA'+DMFT (solid lines) spectra with XPS and BIS ex- 
perimental data (circles, diamonds, crosses) for MnO (upper 
panel), CoO (middle panel) and NiO (lower panel). Fermi 
level is zero. 



high energy structures as Ni-4s and Ni-4p states contri- 
butions. 



Strontium vanadate SrV03 is perhaps one of the most 
simple paramagnetic strongly correlated metallic sys- 
tems. There is no wonder that it is widely used as a test 
system for various LDA+DMFT based numerical tech- 
niques (38l - l4l| . SrVC>3 has ideal cubic perovskite struc- 
ture with one d-electron in V-3d shell within triply de- 
generated t2 9 subshell. LDA and LDA' band structure 
calculations are performed as described in Refs. [38l - l4lT | 
via LMTO method with von Barth-Hedin exchange cor- 
relation energy Q . 

The 3d bands of vanadium cross the Fermi level, while 
oxygen 2p states are at -8 - -2 eV i.e. much lower than 
the Fermi level (see Fig.[7j left panel, dashed lines). If we 
exclude E^® A contribution for V-3d states as described 
in Sec. Ill, we obtain the LDA' band structure shown in 
Fig. [Jj (left panel, solid lines). Similar to Rcf. [l2| within 
LDA' approach energy splitting between V-3d and 0-2p 
bands \Ed — E p \ becomes larger, than in conventional 
LDA. Since the total number of electrons is fixed, the in- 
crease of | Ed — E p | LDA' is related to 0-2p bands going 
down in energy by about 0.5 eV, with V-3d states remain- 
ing almost unchanged. One should mention here also, 
that the overall bandshapes are practically not changed 
in comparison with the conventional LDA bands. The 
same is true of course for densities of states presented on 
the left panel of Fig. [7J 

Another example of paramagnetic strongly corre- 
lated metallic system widely treated by LDA+DMFT is 
Sr 2 Ru0 4 with Ru-4d 4 t 2g subshell (see Ref. gl and ref- 
erences therein). Si - 2Ru04 is a layered perovskite with 
an ideal body-centered tetragonal crystal structure. For 
LDA and LDA' calculations we used settings described in 
Ref. LDA (dashed lines) and LDA' (solid lines) band 
dispersions and DOS'es are plotted in Fig. [7] (right panel). 
The picture here is not that simple as for SrV03. The Ru- 
4d states, crossing the Fermi level, almost preserve their 
energy positions and dispersions within LDA'. However 
LDA' leads to \Ed — E p \ splitting, because of non- uniform 
narrowing of 01-2p and 02-2p states, together with the 
slight shift of 02-2p states. In total \Ed — E p \ energy 
splitting is about 0.5 eV larger for LDA' than in conven- 
tional LDA. 
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FIG. 8: Densities of states calculated with LDA+DMFT 
(dashed lines) and LDA'+DMFT (solid lines) for SrV0 3 : V- 
3d states - upper panel, 0-2p states - lower panel. Fermi level 
is zero. 
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B. LDA+DMFT and LDA'+DMFT DOS 



In contrast to previous works (Refs. |38l - l42j ) we used 
here the full TB-LMTO-ASA calculated LDA and LDA' 
Hamiltonians, employing none of the widely used projec- 
tion techniques. In QMC calculations inverse tempera- 
ture was taken to be j3 — 10cV _1 , with 80 time slices for 
SrV0 3 , while for Sr 2 Ru0 4 we used /3 = 15eV-\ with 64 
time slices. Coulo mb p arameters were taken to be [7=6.0 
eV and J=0.7 eV Hi for SrV0 3 and 3.2 eV and 0.7 eV 
for Sr 2 Ru04 respectively 42j. Number of Monte Carlo 
sweeps was of the order of 10 6 . To obtain DMFT(QMC) 
[29j densities of states at real energies, we again employed 



FIG. 9: Densities of states calculated with LDA+DMFT 
(dashed lines) and LDA'+DMFT (solid lines) for Sr 2 Ru0 4 : 
Ru-4d states - upper panel, 0-2p states - middle and lower 
panels. Fermi level is zero. 



the maximum entropy method [301 ] . To get correspond- 
ing DMFT 0-2p densities of states the method of Pade 
approximants was applied to do analytic continuation 
for DMFT self-energy from Matsubara to real frequen- 
cies, with further crosschecking of "MEM" and "Pade" 
DOS'es to ensure the quality of restored self-energy for 
real frequencies. 

In Figs. [5] and |H] we present the total and par- 
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tial densities of states for SrVC>3 and Sr2Ru04 calcu- 
lated by the conventional LDA+DMFT (dashed lines) 
and LDA'+DMFT (solid lines). For both systems 
LDA'+DMFT results show lower positions of 0-2p states 
in comparison with LDA+DMFT. However, for Sr2Ru04 
it does not reduce just to a rigid shift of oxygen states by 
about 0.5 eV, as in the case of SrVC>3, but is the combi- 
nation of some small shift with non-uniform narrowing of 
oxygen bands. Thus, for Sr2Ru04 only the high energy 
threshold of 0-2p states moves down by 0.5 eV. 

As opposed to Refs. f38U4l| in both calculations for 
SrV03 we observe very smooth upper and lower Hubbard 
bands in V-3d DOS (upper panel of Fig. [HJ . This a gree s 
well with the full orbital calculations reported in Ref. [2J] . 
Also in Ref. [24| it is shown that smaller value of Ed c (if 
Edc is treated as free parameter) moves oxygen states 
down in energy, which leads to better agreement with 
experiment (see the next paragraph). 
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FIG. 10: (Colour online) Comparison of LDA+DMFT 
(dashed lines) and LDA'+DMFT (solid lines) calculated spec- 
tra for SrVOa with experimental data. With cyan colour 
LDA'+DMFT t 2g (solid line) and e g (dot-dash line) V-3d 
contributions are shown. Fermi level is zero. 



C. Comparison of LDA+DMFT and LDA'+DMFT 
results with X-ray experiments 

In Figs. M and CCD LDA+DMFT (dashed lines) and 
LDA'+DMFT (solid lines) calculated spectra for SrV0 3 
and Sr2RuC>4 correspondingly are drawn. To get theo- 
retical spectra from total DOS Gaussian broadening to 
simulate experimental resolution and Lorentzian broad- 
ening to simulate lifetime effects together with multipli- 
cation with Fermi distribution function were performed 
as described elsewhere [38l - l42j . On the figures emission 
(left side) and absorption (right side) spectra are plotted. 

For both systems we have reasonable agreement with 
experimental data (circles) for valence and conducting 
bands HUES! (see Figs. CLS HQ). However strength 
of quasiparticle peak is a bit overestimated for va- 
lence band and underestimated for conduction band in 
both LDA+DMFT and LDA'+DMFT methods. The 
LDA'+DMFT results give slightly better energy position 
of 0-2p states in comparison to LDA+DMFT. In general 
obtained by LDA'+DMFT results are in agreement with 
previous LDA+DMFT works (see Refs. jaiB3|). 

To demonstrate presence of well known lower Hubbard 
band at -1.5 eV for SrV0 3 (a§E| on left panel of Fig. CDS 
V-3d t2 g contribution is shown by cyan line. In Fig. [10] 
(right panel) for SrV03 instead of upper Hubbard band 
around 2.5 eV LDA'+DMFT shows rather broad shoul- 
der. This shoulder is formed by ti g (solid cyan line) and 
e g (dot-dash cyan line) V-3d contributions which corre- 
sponds to previous works [38r44ll |. However the e g sub- 
band in our case is also modified by correlations. It is 
shifted up on about 1 eV (as should be for completely 
empty states) and it has smaller width compared to the 
LDA one. For Sr2Ru04 it is known that correlations 
lead to formation of lower Hubbard band satellite near 
-3 eV H3). This satellite is also seen in the LDA'+DMFT 
results on the right panel of Fig. [TT] and is formed essen- 
tially by Ru-4d t^g states (cyan line). 




Energy, eV Energy, eV 



FIG. 11: (Colour online) Comparison of LDA+DMFT 
(dashed lines) and LDA'+DMFT (solid lines) calculated spec- 
tra for Sr2Ru04 with experimental data. With cyan colour 
LDA'+DMFT t 2g Ru-4d contribution is shown. Fermi level 
is zero. 



VI. 5. CONCLUSION 

This work continues our research of the double count- 
ing problem arising within the LDA+DMFT compu- 
tational scheme. The problem appears because some 
portion of local electron-electron interaction is already 
present in LDA calculations. Since DMFT gives ex- 
act local solution of the Hubbard-like model one should 
avoid double counting between LDA and DMFT local 
electronic interactions. Despite 15 years of developing of 
the LDA+DMFT method still there are no unique defini- 
tion of this double counting term. This happens because 
LDA contribution to exchange correlation energy has no 
diagrammatic expression. Several different ad hoc defi- 
nitions which are available now work well only in some 
particular cases, for some particular compounds. Some- 
times one can get even qualitatively wrong LDA+DMFT 
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solution if double counting term is chosen not careful 
enough. To overcome this problem we proposed consis- 
tent LDA'+DMFT approach [l2|. It uses natural as- 
sumption of explicit exclusion of LDA exchange correla- 
tion potential for correlated electronic shells since any- 
how exchange-correlation effects will be accounted later 
by DMFT. Then local interactions left out for correlated 
states in the LDA' are only Hartree ones. After that cor- 
responding double counting term of the LDA'+DMFT 
Hamiltonian consistently must be taken in the local 
Hartree form (FLL form). 

With this paper we present extensive LDA'+DMFT in- 
vestigation of typical representatives of two wide classes 
of strongly correlated systems in the paramagnetic phase: 
strongly correlated metals (SrVC>3 and S^RuCU) and 
charge transfer insulators (MnO, CoO and NiO). For 
strongly correlated metals where double counting is not 
that severe LDA'+DMFT agrees well with traditional 
LDA+DMFT results with FLL double counting type. 
LDA'+DMFT gives slightly better position of 0-2p 
states in comparison with experiment. LDA'+DMFT re- 
sults for charge transfer insulators MnO, CoO and NiO 
are more interesting. CoO and NiO systems arc found 
to be metals within conventional LDA+DMFT calcu- 
lations while LDA'+DMFT gives proper insulating so- 
lution. Transition metal 4s-states missed in previous 



LDA+DMFT works on these monooxides are found to be 
responsible for charge gap asymmetry around the Fermi 
level. 

Finally one can conclude that proposed by us consis- 
tent LDA'+DMFT method works well for both metal- 
lic and insulating systems. We believe that our 
LDA'+DMFT provides reasonable parameter free treat- 
ment of the double counting problem. 
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TABLE I: LDA and LDA' occupancies and corresponding values of LDA+DMFT and LDA'+DMFT double counting terms 
(eV) for systems under consideration. 



Compound 


IlLDA 




LDA+DMFT 
FLL(LDA) 


LDA+DMFT 
FLL(SC) 


LDA'+DMFT 
FLL(LDA) 


LDA'+DMFT 
FLL(SC) 


S1-VO3 


2.61 


2.44 


12.33 


11.99 


10.35 


10.92 


Sr 2 Ru0 4 


5.65 


5.39 


14.32 


14.60 


12.92 


13.73 


MnO 


5.59 


5.43 


39.05 


35.49 


36.62 


35.30 


CoO 


7.60 


7.41 


54.28 


50.90 


51.42 


50.49 


NiO 


8.54 


8.34 


60.90 


62.01 


57.91 


58.13 



